Negative differential resistance due to the resonance coupling of a quantum-dot dimer 
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Electron tunneling through a coupled quantum-dot dimer under a dc-bias is investigated. We find 
that a peak in the I-V curve appears at low temperature when two discrete electronic states in the 
two quantum dots are aligned with each other - resonance coupling. This leads to a negative differ- 
ential resistance. The peak height and width depend on the dot-dot coupling. At high temperature, 
the peak disappears due to thermal smearing effects. 

PACS numbers: 73.23.-b, 73.63.Kv 
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I. INTRODUCTION 

A large number of studies have focused on nanostruc- 
ture materials, such as quantum dots (QDs) 1], as we 
move into the era of nanoscience and nanotechnology. 
This is due largely to academic interests and their po- 
tential applications. It has been proposed that quantum 
dots are used as building blocks of electric circuits^ and 
even of quantum logic gates 0. Since most of these build- 
ing blocks involve many QDs, many experimental0,0,0 
and theoretical0,H,H investigations are focused on more 
than one QD system, especially two-QD system - QD 
dimer. In a two QD system, not only the charging ef- 
fect but also the interdot coupling and the alignment of 
electronic states in the two QDs play important roles. 
Because of the interdot coupling, electrons can be shared 
by both QDs. A state similar to the covalent state in a 
molecule can be formed in a QD dimer. This state can 
be manipulated by varying the external parameters of 
the system, such as the interdot coupling[6|]. Therefore 
QD dimers are proposed to be candidates for building 
quantum logic gatesQ- The interdot coupling also yields 
new features in the Coulomb Blockade conductance spec- 
troscopy. For example, a system with two isolated iden- 
tical QDs has an oscillation structure in the conductance 
spectroscopy in the Coulomb blockade regime. By tun- 
ing the interdot coupling, those peaks are split into two 
peaks each.0, 0. The alignment of electronic states in 
two QDs makes the I-V characteristics of a dimer with 
two QDs in series much different from that of a single 
QD system with a step-like structure. The I-V curve of 
a QD dimer has many peaks and these peaks are due to 
the resonant tunneling when two electronic states in the 
two QDs are aligned[a,|3- 

Most of the previous studies are on a QD dimer with 
two QDs in series. In this paper, we shall study a system 
with two QDs coupled in parallel with source and drain 
leads. One QD (QD1) is connected to both source and 
drain leads while the other QD (QD2) is connected to 
only one lead. This configuration may arise when one 
performs STM experiments on quantum dots on a sub- 
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strate. We show that the alignment of electronic states 
in the coupled QD dimer, which we shall call resonance 
coupling, can modify greatly the electron transport. We 
find that a peak in I-V curve occurs when two electronic 
states in the two QDs are aligned, leading to a nega- 
tive differential resistance (NDR) at low temperature. In 
the presence of the electron-electron (e-e) interaction, the 
peak splits into two. 

The present paper is organized as follows. Our model 
and the method used to calculate the I-V curves are de- 
scribed in Sec. |nj Then, we present, in Sec. II I II our 
calculated results for different cases. The numerical ev- 
idence of a NDR in our model is given. We shall also 
provide the explanation of the occurrence of the NDR. 
The summary in Sec. IIVI concludes the presentation. 



II. MODEL AND METHOD 

As shown in Fig. ^ we consider a QD dimer with two 
coupled QDs connected to metallic leads in parallel. One 
of the leads act as a source, say the left lead (L-lead) in 
the figure, and the other (the right lead) as a drain. One 
quantum dot (QD1) is connected to both external leads 
while the other quantum dot (QD2) is only connected to 
one of these leads, say the right lead (R-lead) as shown 
in the figure. Thus, electrons, flowing from the L-lead to 
the R-lead through QD2 must tunnel to QD1 first. This 
set-up can describe an STM tunneling experiment of a 
QD (QD1) on a substrate while there is an adjacent QD 
(QD2) coupled to QD1. 

We shall model both the L- and R-leads as one- 
dimensional (ID) non-interacting Fermi gases. We as- 
sume that the electron level spacing in both QDs is large 
in comparison with the level broadening due to the tun- 
neling process and the thermal effect. Thus, the spectra 
of both QDs are discrete. For simplicity, we consider 
that there is only one energy level with spin degeneracy 
in each QD. The electron-electron (e-e) interaction in a 
mesoscopic system is usually larger than or comparable 
to other energy scales like the typical level spacing and 
the thermal energy. Hence the Coulomb blockade effect 
plays an important role in electron transport when QDs 
are weakly coupled to external leads. We shall consider 
the e-e interaction in our model. And we shall describe 
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FIG. 1: Schematic diagram of a system with two coupled 
quantum dots connected to source and drain leads in parallel. 
One quantum dot (QD1) is connected to both the left (L) and 
the right (R) external metallic leads. The other quantum dot 
(QD2), is only connected to R-lead. 

the tunneling process between leads and QDs with a tun- 
neling Hamiltonian. 

The total Hamiltonian of this system can be expressed 

as 

H = H L , R + H C + H T , (1) 

where Hl,r is the Hamiltonian of two ID ideal metallic 
leads, H c is the Hamiltonian of the central region with 
two coupled QDs, and Ht is the tunneling Hamiltonian. 
The Hamiltonian of two ID ideal metallic leads is 

Hl.r = Y %kcl ka c lM , (2) 

where a =|, I is the spin index, and . (c t fc ) is the 
creation (annihilation) operator of an electron with en- 
ergy e^k and spin a in lead I = L, R, the left or right 
leads. 

The Hamiltonian of the central region is 

H c = Y +Y( td ^ + H - c -) 

i=l,2 cr 

+ Y Uin^n^, (3) 

i=l,2 

where d\ a (d i(T ) is the creation (annihilation) operator of 
an electron with energy ej and spin a in QD i. t is the 
hopping energy between the states ei and 62 in the two 
QDs. The terms containing Ui describe the e-e interac- 
tions in QD i. 

The tunneling Hamiltonian is given as 

H T = Y V lM , i 4 ika d ia +U.c, (4) 

l , fc ,CT 

i = l,2 

where Vi ka,i is the coupling constant between state e;^ 
in the lead I and state €i in the QD i. Here we assume 
that the tunneling process conserves electron spins. 



By using the non-equilibrium Green function method, 
the current for a steady state is given as[ToL ITlj. 

/ = E^/ dETr(^<G>-^>G<), (5) 

where G„ is the lesser (greater) Green function of the 
central region with spin a. Y%' < ^ > ' is the lesser (greater) 
self-energy matrix of the spin a electron in the central re- 
gion with the contribution from L-lead only. Below, we 
shall briefly outline the steps of calculating these Green 
functions. We shall see that the self-energies can be 
treated as parameters and their values can be chosen ac- 
cording to our considerations. 

The lesser (greater) Green function can be calculated 
by the Keldysh equations, 

G<<» = (6) 

where G£ is the retarded (advanced) Green function 

of the central region with spin a. E^ is the lesser 
(greater) self-energy of spin a electron in the central re- 
gion with the contribution from both leads. 

The matrix elements of the retarded and advanced 
Green functions are defined as 

Gl tij (t-t') = -iO(t-t')({d ia (t),dl(t')}), (7a) 
Gl M {t - f) = i6(t' - t)({d i(T (t),dl(t')}), (7b) 

where 0(t) is the step function, and {a, b} — ab + ba is 
the Fermion anti-commutator. 

The equation-of-motion method is used to calculate 
the retarded (advanced) Green function, G£ . The 
time derivative of the retarded (advanced) Green func- 
tion Ga^ is 

■^ G r(a) = gQ. _ + AG r(a) + UG r(u)(2) + jga) ^ (g) 

where 

and / is the unity matrix. E^ a ' 1 is the retarded (ad- 
vanced) self-energy with spin a. The elements of the 
second order Green function G^ " 2 ' are defined as 

<? 2) = Tie(±tTt'){{d i(T (t)n m (t),dl(t')}), (10) 

where a is the opposite spin to er, and n is = d\ s d is . 

Usually the equation-of-motion of the Green func- 
tion Ga^ contains higher order correlation functions, 
for example, Gl {aK2} in Eq. ©. Whc n we calculate 
the equations-of-motion of these higher order correla- 
tion functions, even higher order correlation functions 
appear. In order to close those equations, we must trun- 
cated those equations at an appropriate order. Here, we 
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only keep those equations containing terms of correlation 
functions up to the second order. Then, we can get the 
expression of G^^ by solving these equations. Since the 
calculation process is well established, we don't give it 
explicitly here. The detail calculation can be found in 
Ref. [13 ■ 

The retarded (advanced) Green function Ga^ con- 
tains terms with the average electron numbers of the 
opposite spin, (n is ) (i — 1,2), in QD i. The average 
electron number in QD i can be calculated by 

/dE 
— lmG< u (E). (11) 

Thus, we need to calculate them self-consistently. After 
evaluating the average electron numbers in both QDs self- 
consistently, we can obtain all Green functions. Because 
we consider here a system with two QDs, all these Green 
functions and self-energies are 2x2 matrices. 

As shown above, we need to know self-energies to cal- 
culate those Green functions. The elements of lesser 
(greater) self-energy are given as 

v<<» v* „*><(>) 

l,k 

= S^ (>) + S%< (>) , (12) 

where 

S ^ (>) = E v^.iVtVX* <(>) > ( 13a ) 

k 

^« (>) =E^.*V3,^9& <(>) , (13b) 

k 

are the lesser and greater self-energies with contributions 

from both L- and R- leads, respectively. g& is lesser 
(greater) Green function for lead I with spin a and the 
elements of these two Green functions are given as 

= ifM,k)e- Ul ^ t - t '\ (14a) 

&-t') = ^(c Lk Jt)cl ka (t>)) 

= -i(l-fl(ei, k ))e- Ul ^ t - t '\ (14b) 

where /z(e) = (1 + e /3 ' e_M ^) _1 is the Fermi-Dirac distri- 
bution function, fii is the chemical potential in lead I 
and (3 — l/(ksT). The time dependent operators in lead 
I are c lM (t) = e-^\ ka , c\ M (t) = e^c\ M . 

The elements of the retarded (advanced) self-energy 
are given as 

l,k 

where g l J r ^ is the retarded (advanced) Green function 
of lead I with spin a and their elements are given as 

9 l ;f\t ~ = T*0(±t T t')({c lM (t), c\ M {t')}) 

= T i6(±tTt')e-' lei '« < - t - t '\ (16) 



In our calculation, we need to know the Fourier trans- 
formation of the elements of the retarded (advanced) self- 
energy. They are 

l,k 

= K%](E) T ^(E). (17) 

where the real and imaginary parts, that is, the level-shift 
function Kl {a \E) = A^ r(a) {E) + K RM<1 \e), and the 
level-width function T r a {a) {E) = T L a Ha) {E) + T RAa) (E) 
are due to the tunneling process between the QDs and 
both L- and R-leads. Usually, these functions are energy- 
dependent. However, under the wide band approxima- 
tion, they do not depend on electron energy. We shall use 
this wide band approximation and assume the level-shift 
and level-width functions to be independent of energy. 
In the absence of a magnetic field, they do not depend 
on spin index cr, and we shall drop it in our notations 
below. Therefore, we can define those self-energies as 
parameters with four 2x2 constant matrices A L ( R *> and 
Y L ( R ) instead of using the coupling constant {Vi t ka,i}- 
Then the retarded (advanced) self-energy is 

According to Eq. i|12|) . the lesser self-energy is given as 

E<(E)=i(T L f L (E)+T R f R (E)) (19) 

and the greater self-energy is 

£>(25) - -*(r L (l - f L (E)) + V R {1 f R (E))). (20) 

The lesser and greater self-energies depend on en- 
ergy only through the Fermi-Dirac distribution function 

fh{R){E)- 

The Green functions can be calculated by equation-of- 
motion method described above in terms of all the self- 
energies. Then we can evaluate the I-V curve by using 
Eq. J5J . In the next section, we shall give some numerical 
results of I-V curves for different cases. 



III. NUMERICAL RESULTS AND DISCUSSION 

A. In the absence of e-e interactions in both QDs 

First, we consider a simple case that the electron- 
electron (e-e) interactions in both QDs are absent, i.e., 
U\ = U2 = 0. The energy needed to add an extra elec- 
tron in a QD due to the e-e interaction is e 2 /C, where 
C is the capacitance of the QD. When the size of a QD 
are large, this energy may be very small compared to the 
level spacing. We can then ignore the e-e interactions in 
the QD. Thus, this may be used to describe a relative 
large QD. 
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In order to investigate the resonance coupling effect on 
the I-V characteristics, we assume that as shown in Fig.Q] 
one of two QDs (QD1) is connected to both external leads 
while the other one (QD2) is only connected to the right 
lead. When no external bias is applied, electronic state 
energy ei(0) in QD1 is set to be smaller than energy £2(0) 
in QD2. The external bias is applied by increasing chem- 
ical potential /xl of the left lead and keeping fj,R of the 
right lead unchanged. We assume that the barriers be- 
tween QD1 and left and right leads are symmetric, and 
voltage drops uniformly across the QD1. Thus, the en- 
ergy level ei of QD1 shifts to t\ — ei(0) + eV/2 under a 
bias V bias. In reality, because of the difference of the 
electro-static potentials of QD1 and R-lead, there is a 
voltage drop in QD2 and energy £2 in QD2 should also 
shift. However, t\ shifts more than £2 does. When a cer- 
tain bias is applied, £1 and £2 can be tuned to be aligned 
with each other, or resonance coupling. In fact, other 
model parameters may also change. It is known that bias 
dependence of other model parameters may also affect I- 
V characteristics However, for simplicity, we assume 
that the energy level £2 in QD2 remains unchanged be- 
cause we are interested in the resonance coupling effect 
on the I-V curves in this study. This assumption may 
affect the positions of / — V peaks mentioned below, it 
shall not change the physics studied. We shall come back 
to this point in our discussions. 

Accordingly, we set the energies of two electronic states 
in the two QDs at zero bias, £i(0) = 1.0, £2(0) = 3.5. 
Since QD1 are connected to both leads and QD2 are only 
connected to the right leads, we set rf\ = r{\ = 0.2, 

r£ 2 = 0.0, rf 2 = 0.2, rf 2 = r& = o.o, rf 2 = r& = 

0.1. Because the level-shift functions A L ^ only add a 
constant energy shifting to the electronic state in each 
QD, their values don't change our final results. We set 
them to zero, that is, A^- = 0.0, where i, j — 1, 2. The 
chemical potentials of the left lead fiL — eV can vary 
while /if/ = 0. All these parameters are in an arbitrary 
unit of energy T which is a system parameter. 

We first calculate I-V curves at a very low tempera- 
ture, T = 0.01 in the unit of Y/ks, so that we can neglect 
the thermal broadening of energy levels in QDs. Fig.^a) 
shows the calculated I-V curve and Fig. |5{b) shows the 
numerical result of the bias dependence of the average 
electron numbers (n^ in QD i. At low and high bias, 
the two electronic states are far apart. Since QD2 is con- 
nected to the right lead only, electrons can only tunnel 
into it through QD1. In the absence of any inelastic pro- 
cesses, the chance for an electron tunneling from QD1 to 
QD2 is very small when the two electronic states, t\ and 
£2, are far apart. Therefore, the average electron number 
in QD2 is almost zero. Its effect on the tunneling process 
is negligible, and the I-V characteristics should be very 
similar to that of a single QD system. This is indeed 
what we can see from Fig. Ela). As shown in Fig. Ufa), 
there is a peak at V — 5.0, where two electronic states 
£1 and £2 are aligned with each other. Correspondingly, 
the average electron number in QD1, (n^) has a valley at 



this bias as shown in Fig.JSfb). At the same time, there 
is a dramatical increase in the average electron number 
in QD2. This is not surprising at all. The chance for an 
electron tunneling from one electronic state in QD1 to 
another state in QD2 increases as the energies of the two 
states approach each other. It become maximum when 
they are precisely aligned. We shall call it resonance cou- 
pling of the two QDs when this happens. At resonance 
coupling, a new tunneling channel through QD2 opens. 
In another words, the effective tunneling rate of an elec- 
tron out of the QD dimer increases. This causes the peak 
observed in Fig.^a), leading to a NDR. 
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FIG. 2: Numerical results of I-V and {n^-V in the absence of 
e-e interactions. The parameters are ei(0) = 1.0, £2(0) = 3.5, 
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0.2, li 2 = 0.0, lf 2 = 0.2, rf 2 = Iai = 0.0, 



rf a = rfi = 0.1, Af} R) = 0.0 = 1,2), n t = eV and 
fiR = 0, where all these parameters are in unit of F. The 
temperature T — 0.01 in unit of F/fcs. (a) Current I vs 
applied bias V. The current I has a peak at V = 5.0. (b) 
Average electron numbers in both QDs, (n 4 ) vs applied bias 
V. (wi) has a valley while (n 2 ) has a peak at V = 5.0. 



This can be understood from the following argument. 
As the electron tunneling probability from QD1 to QD2 
increases, the effective tunneling rate out of the electronic 
state £1 increases and the effective tunneling rate into the 
electronic state £2 increases as well. Therefore, (n x ) de- 
creases and (n 2 ) increases because more electrons tunnel 
from QD1 to QD2. For the sequential tunneling, the res- 
onant tunneling current and the average electron number 
in the resonant level at zero temperature are given as|llj 
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where r in ' out ) is the tunneling rate of an electron into 
(out of) the resonant level. Near the resonance coupling, 
for QD1, the tunneling rate from the left lead into 
the dot remains almost unchanged and the tunneling rate 
pout 0lr (. Q f dot increases because of a new tunneling 
channel through QD2. For QD2, the tunneling rate T 2 n 
from QD1 into QD2 increases while the tunneling rate 
r 2 ut out of QD2 into the right lead remains unchanged. 
Therefore, the overall tunneling rates for the QD dimer is 
as follows. The tunneling rate r in from left lead into the 
QD dimer is almost unchanged, while the effective tun- 
neling rate r out out of the QD dimer increases. There- 
fore, the total current increases and (n^ decreases while 
(n 2 ) increases according to Eq. (|21|l . In turn, it generates 
a I-V peak and a NDR. 
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FIG. 3: Numerical results of I-V curves at different hopping 
energies and temperatures. All other parameters are the same 
as those in Fig. [5] (a) I-V at different hopping energies, 
t = -0.05 (solid line), t = -0.10 (dashed line), and t = -0.30 
(dotted line). The peak height and width increase as the 
hopping energy t. (b) I-V at different temperatures, T — 0.01 
(solid line), T = 0.50 (dashed line), and T = 1.00 (dotted 
line). The peak is smeared at high temperatures. 

We have seen that resonance coupling between two 
dots can lead to a peak in the I-V curve. Thus, we 
should expect the height and width of this peak to de- 
pend on the interdot coupling strength. Also, the ther- 
mal energy can introduce inelastic tunneling processes 



which can wash out the resonance effect. One should 
also expect that the peak is also sensitive to the tem- 
perature. The I-V curves at different interdot coupling- 
strengths between two QDs, and at different tempera- 
tures are shown in Fig. O^a) and (b), respectively. As 
shown in Fig. |3 the peak height and width increase as 
the coupling between two QDs increases. Because the 
chance of an electron tunneling from QD1 to QD2 in- 
creases with the interdot coupling, more electrons can 
tunnel out of the dimer through QD2, and the effective 
tunneling rate out of the QD dimer increases. There- 
fore the peak height and width increase as the interdot 
coupling increases. We can also see the thermal smear- 
ing of the NDR. The peak disappears gradually with the 
increase of temperature. 



B. The presence of e-e interactions in QDs 

In order to see whether the I-V peaks, thus the NDR, 
due to the resonance coupling of two QDs will survive 
when the e-e interaction is present, we repeat the above 
calculation by including non-zero U\ and [/ 2 . We set 
the Coulomb interactions Ui = U2 = 5.0 while keeping 
the other parameters the same as those in Fig. [21 The 
calculated results of I-V curve and the bias dependence 
of the average electron numbers (n,) in QD i are shown 
in Fig.^Ja) and (b), respectively. As one will expect, the 
peak at V = 5.0 is still there with the similar features as 
those in the absence of the e-e interaction. The reason 
for the occurrence of this peak should be the same as 
that in the previous section. 

However, there is an extra peak at V = 15.0. This new 
peak can be attributed to the e-e interactions in QDs. In 
order to understand the origin of this new peak, let us 
examine the possible electronic states in QD i. In the 
case that no electron in the QD, an electron of energy 
e; can move in freely into it. However, only the electron 
with energy ei + Ui can jump into the QD when there 
is already one electron in because of the e-e interaction. 
Therefore, the electronic states of the QD has effectively 
two distinct energy values, Ci and Ci + Ui. We have al- 
ready identify the first peak as the resonance coupling 
due to the alignment e\ and e 2 . It is natural to expect 
another peak to appear when e\ is aligned with e 2 + J7 2 - 
In our case, it corresponds to V — 15.0, exactly what 
we see in Fig. EJa). In comparison with the first peak, 
the second peak is smaller and narrower. This is due 
to the different properties of the two electronic states e 2 
and e 2 + Ui- Unlike e 2 , e 2 + U2 doesn't exist unless one 
electron has already been in QD2. However, when t\ and 
(-2 + U2 are aligned, e 2 and ei are far apart. Therefore, 
the probability of an electron hopping from ei to e 2 is 
small. Consequently, the average electron number in e 2 
is very small. Thus, the existing probability of e 2 + U2 
state is very small too. This is probably the reason why 
the second peak is smaller and narrower. 

Unlike the first resonance-coupling peak where (n 2 ) is 



6 




12 16 20 



0.3 

u 
u 

I 

Z0.2 
c 
o 
u 

o 
u 

■20.1 



< 0.0 



V(r/e) 



(b) 










\ 

<nj> 






<n 2 > 



12 16 20 



V(r/e) 



at the second peak. 

In the above calculations, we have not considered the 
case that t\ + U\ is aligned with e 2 + t/ 2 . This is be- 
cause the bias is too small to allow electrons tunneling 
through ei + U\, when those two electronic states are 
aligned. Moreover, if we set the same e-e interactions in 
two QDs, t\ + U\ and e 2 + U2 are aligned at the same 
time when e± and e 2 are aligned. These two effects may 
not be distinguished easily. In order to avoid this pos- 
sible confusion, we consider different e-e interactions in 
two QDs such that e\ and e 2 are not aligned with each 
other when ei + U\ and e 2 + U2 do. We also choose a 
set of parameters in such a way that electrons can tunnel 
through ei + U\ when it is aligned with e 2 + J7 2 . 




FIG. 4: Numerical results of I-V and (n^-V in the case in 
presence of e-e interactions. The e-e interactions in two QDs 
are Ui = U2 = 5.0 and other parameters are the same as those 
in Fig. [5] (a) Current / vs applied bias V. The current I has 
two peaks, one is at V — 5.0 and the other is at V — 15.0. 
(b) Average electron numbers in both QDs, (n^) vs applied 
bias V. (%) has a valley at V = 5.0 and a peak at V — 15.0. 
(712) has two peaks. One is at V — 5.0 and the other is at 
V = 15.0. 



maximum and (n^ is minimum, both (n^ and (n 2 ) are 
maximum at the second peak. (n 2 ) becomes maximum 
at the second resonance-coupling peak because the effec- 
tive tunneling rate into QD2 from QDI increases as ei 
approach e 2 + J7 2 . The maximum of (n x ) is because both 
tunneling rate into and out of QDI increase. Around 
V = 15.0, the effective tunneling rate out of QDI in- 
creases because more electrons can tunnel from QDI to 
QD2. On the other hand, the effective tunneling rate 
into QDI increases too. Two electrons can now be in 
QDI at the same time, occupying ei and e± + U\ states, 
respectively. After the electron with energy e\ tunnels 
out of QDI, the remaining electron has energy t\ instead 
of its original energy e\ + U\. Then a new electron with 
energy t\ + U% can tunnel from the left lead into QDI. 
As e% and e 2 + U2 are aligned, electrons in e\ tunnel out 
of QDI faster. Therefore, electrons can tunnel into QDI 
from the left lead faster and the effective tunneling rate 
into QDI increases. So when e\ and e 2 + 1/ 2 are aligned, 
both the effective tunneling rate into QDI P" and the 
effective tunneling rate out of QDI T° ut increase. Ac- 
cording to Eq. 121|l the average electron number in ei 
may increase if T™ increases faster than T° ut . This ex- 
plains the possibility that {n^ reaches maximum value 
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FIG. 5: Numerical results of I-V and {n^-V in the case of 
asymmetric e-e interactions of two QDs. The e-e interactions 
in two QDs are Ui = 3.0 and U2 ~ 5.0. All other parameters 
are the same as those in Fig. [5] (a) Current / vs applied bias 
V. The current / has four peaks, one is at V = 5.0, one at 
V = 15.0, and other two small peaks at V = 8.3 and V = 8.6, 
respectively, (b) Average electron numbers in both QDs, (rij) 
vs applied bias V. (ji-i) has a valley at V = 5.0 and three 
peaks. The first is at V = 8.3, the second at V = 8.6, and the 
third at V = 15.0. (n 2 ) has three peaks. One is at V — 5.0, 
one at V = 8.6, and the last at V — 15.0. 

Fig. |S{a) shows I-V curve with U\ = 3.0, U2 = 5.0, 
and all other parameters are the same as those in Fig. [21 
Besides two peaks at V = 5.0 and V = 15.0, there are 
two additional very small peaks. One is at V = 8.3 and 
the other at V = 8.6. Similar to the two old resonance- 
coupling peaks, these two new small peaks are due to the 
alignment of ei + U\ and e 2 + C/ 2 . This resonance-coupling 
yields two peaks instead of one because of the competi- 
tion of two opposite effects that do not exist in the pre- 
vious case. As we explain above, the existing probability 
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of electronic state ti + U2 decrease as e\ and 62 are far 
apart. However, when ei + U± and £2 + U2 become closer, 
the electron tunneling probability from e% + U± to e 2 + Ui 
increases. The competition of these two effects leads to 
those two peaks in the I-V curve. 

As shown in PigEfb), (%) has one valley at V = 5.0 
and three peaks. They are located at V = 8.3, 8.5, and 
15.0, respectively. (n 2 ) has three peaks at V — 5.0, 8.5, 
and 15.0, respectively. The behavior of {nj, i = 1,2, at 
V = 5.0 and 15.0 has been already explained before. The 
reason of the occurrence of the other two peaks of (%) is 
due to the competition of two opposite effects mentioned 
above. When ei + U\ and e 2 + U% are aligned, there 
are those two opposite effects. One decreases the average 
electron number in QD2 while the other increases it. But 
the increase of the average electron number is larger than 
the decrease of it, therefore we can only see one peak of 
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FIG. 6: Numerical results of I-V curves at different coupling 
strengths and e-e interactions. All other parameters are the 
same as those in Fig. El ( a ) I-V at different hopping energies, 
t = -0.05 (solid line), t = -0.10 (dashed line), and t = -0.30 
(dotted line). The e-e interactions are U\ — U2 — 5.0. (b) 
I-V at different e-e interactions, Ui = U2 — 3.0 (solid line), 
Ui = U 2 = 5.0 (dashed line), and Ui = U 2 = 1.00 (dotted 
line). 

The I-V curves of different coupling strengths and e-e 
interactions are shown in Fig. Ela) and (b), respectively. 
As one will expect, the heights and widths of both peaks 
increase as the coupling strength between two QDs in- 
creases. From Fig. Efb), one can see that the e-e inter- 
actions have no effect on the peak at low bias because 
the peak at low bias is due to the alignment of ei and 
£2- A strong e-e interaction shifts the position of the 



peak at high bias, but it does not change its height and 
width. The peak at high bias is due to the alignment of 
ei and £2 + C^2- As we have explained before, the height 
and width of this peak depend on the existing probabil- 
ity of £2 + U2 state. When ei and £2 + U2 are aligned, ei 
and 62 are far apart. Thus the probability of an electron 
hopping from ei to €2 is too small to be sensitive to the 
change of the e-e interactions. Consequently, the exist- 
ing probability of £2 + U2 state is independent of the e-e 
interactions. Therefore, the strong e-e interactions have 
little effect on the height and width of this peak. The 
I-V curves at different temperatures are shown in Fig.0 
Peaks disappear at hight temperature due to the thermal 
smearing effects. 
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FIG. 7: Numerical results of I-V curves at different temper- 
atures, T = 0.01 (solid line), T = 0.50 (dashed line), and 
T = 1.00 (dotted line). All other parameters are the same as 
those in Fig. Peaks are smeared at high temperatures. 

What is seen is a new mechanism for the NDR in a 
system with more than one QDs. The NDR is a very 
important phenomenon and has many important appli- 
cations. Devices with the NDR have been widely used to 
make amplifier and oscillators in a very wide frequency 
range [liij. In superlattice, it is known that a NDR leads 
to many interesting phenomena, such as current-voltage 
oscillation on the sequential resonant tunneling plateau, 
current self-oscillation, and chaos^^. There are many 
mechanisms for NDR, such as Gunn effect and resonant 
tunneling in superlattices. The NDR found here is due to 
the resonance coupling between two QDs. As shown in 
Fig. 03 the wave function of the bonding state is largely 
localized in QD1 in off-resonance case. Electrons can 
only tunnel from the source to the drain through QD1. 
In resonance coupling case, electrons in the bonding state 
can be in both QDs, and electrons can also tunnel from 
the source to the drain through QD2. A new tunneling 
channel is open and a peak can be observed in the I-V 
curve. This leads to the NDR. Furthermore, the proper- 
ties of the NDR depend on the properties of both QDs. 
For example, the widths and heights of peaks vary as the 
dot-dot coupling. And the positions of peaks depend on 
the energy spectra of both QDs. The variations of the 
energy spectra of the QD under the probing tip or the 
adjacent QD, or both QDs may change the positions of 
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FIG. 8: Schematic diagrams of the wave functions of an elec- 
tron in (a) off-resonance coupling and (b) resonance coupling 
cases. <pi and 952 are eigenfunctions with energy ei and £2 in 
two QDs without interdot coupling, ip is the wave function 
(bonding state) in the presence of the coupling, (a) In off- 
resonance coupling case, the overlap of ip\ and ifi2 in two QDs 
is negligible and an electron in the bonding state is largely 
localized in QD1. (b) In resonance coupling case, the over- 
lap of ip\ and <p2 are very large and electron in the bonding 
state has approximately the same probability of being found 
in either one of the two QDs. 

peaks in I-V curves. Thus, the I-V characteristics shall 
have different behavior when the tip is above QD2 rather 
than QD1. Comparing this difference will allow one to 
distinguish a NDR due to the mechanism proposed here 
from the others. A better understanding of this reso- 
nance coupling effect may also enhance the STM as a 
powerful probe not only for a regular surface, but also 
for a cluster structure. 

We show that an I-V peak, thus NDR, appears at the 
resonance coupling between two QDs where their energy 
levels are aligned. For simplicity, we consider essentially 
only one electronic state in each QD. The mechanism 
should survive in the case that there are more than one 
electronic state in each QD. This effect should be pro- 
nounced when the sizes of QDs are small such that the 
discrete natural of electronic state are clear. We assume 
also that the bias shifts only the energy levels in the 



QD directly connected to two leads and has no effect on 
other parameters. In reality, the bias not only shifts the 
ener gy le vels of both QDs but also changes the tunneling 
rates |l2j. But in our quantum dot configuration, which 
may arise in a STM experiments, the energy levels in 
QD1 shall be affected more than those in QD2 by an 
external bias. Thus, the resonance coupling will occur 
under certain bias. We use this simplified assumption in 
order to unambiguously identify the observed I-V peaks. 
Since the physics is due to the resonance coupling which 
indeed should occur in real experiments, we believe this 
new mechanism is very robust, and does not depend on 
our assumption. Of course, the peak positions, which 
are not our main concern, depend surely on the assump- 
tion. In fact, the NDR due to the resonance coupling 
of two QDs may have been observed already in a recent 
experiment |15| . 

IV. SUMMARY 

In summary, the dot-dot coupling effect on the tunnel- 
ing current under a bias is investigated. We considered a 
QD dimer connected to two metallic leads in such a way 
that one QD in the dimer is connected to the both leads 
while the other is only directly connected to one lead. 
We show numerically that peaks, thus NDRs, appear in 
the I-V curves at resonance coupling when two electronic 
states in two QDs are aligned. This NDR exists both in 
the absence and in the presence of e-e interactions. We 
study how the peaks in the I-V characteristics change 
with interdot coupling, e-e interaction, and the tempera- 
ture. We show that the heights and widths of those peaks 
increase with the interdot coupling of the QD dimer. But 
they are almost unaffected by the e-e interactions in both 
QDs when the e-e interactions are strong. We show that 
those peaks are smeared at high temperatures due to the 
thermal smearing effects. 
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